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Abstract. We study the properties of wave functions and the wave-packet 
dynamics in quasiperiodic tight-binding models in one, two, and three dimensions. 
The atoms in the one-dimensional quasiperiodic chains are coupled by weak and 
strong bonds aligned according to the Fibonacci sequence. The associated d- 
dimensional quasiperiodic tilings are constructed from the direct product of d 
such chains, which yields either the hypercubic tiling or the labyrinth tiling. This 
approach allows us to consider rather large systems numerically. We show that 
the wave functions of the system are multifractal and that their properties can 
be related to the structure of the system in the regime of strong quasiperiodic 
modulation by a renormalization group (RG) approach. We also study the 
dynamics of wave packets to get information about the electronic transport 
properties. In particular, we investigate the scaling behaviour of the return 
probability of the wave packet with time. Applying again the RG approach we 
show that in the regime of strong quasiperiodic modulation the return probability 
is governed by the underlying quasiperiodic structure. Further, we also discuss 
lower bounds for the scaling exponent of the width of the wave packet and propose 
a modified lower bound for the absolute continuous regime. 



PACS numbers: 71.23.Ft, 72.15.-v 



Submitted to: J. Phys.: Condens. Matter 



Wave Functions, Quantum Diffusion in Quasiperiodic Tilings 



2 



1. Introduction 

In the 1970s works by Penrose and Ammann showed that the Euclidean space can 
be filled gapless and non-overlapping by two or more tiles which are arranged in a 
nonperiodic way according to matching rules |T|. These quasiperiodic tilings showed 
a hierarchical structure and long-range oriental order. The study of such tilings was 
strongly intensified after the discovery of quasicrystals by Shechtman et al. in 1982 [2]. 
These materials are characterized by a perfect long range order without having a three- 
dimensional translational periodicity. The former is manifested by the occurrence of 
sharp spots in the diffraction pattern and the latter in the occurrence of rotational 
symmetries forbidden for conventional crystals pUH]. While the diffractions patterns 
of quasicrystals were not explainable by the structure of crystals with their periodic 
repetition of unit cells, the mathematical concept of quasiperiodic tilings was suitable 
to describe the structure of this new material class [5]. 

Various experimental investigations revealed rather exotic physical properties. 
For instance quasicrystalline surfaces are anti-adhesive in combination with a high 
level of hardness making them suitable for the production of coatings of cookware, 
engines, tools, etc. [6 ]. Further, they possess a rather low thermal and electrical 
conductance although they contain a high amount of well-conducting elements [7JIH]- 
This motivated extensive research to obtain a better theoretical understanding of the 
structure and the physical properties of these materials. Today several exact results 
are known for one-dimensional quasicrystals [9"Mll|. However, the characteristics in 
two or three dimensions have been clarified to a much lesser degree and are mainly 
based on numerical studies. Due to the lack of translation symmetry in quasicrystals, 
one has to study very large approximants to obtain an insight into the physical 
properties of macroscopic quasicrystals. In addition many physical properties in 
quasiperiodic systems show a power-law behaviour with respect to time or system 
size, which indicates the absence of a typical time/length scale in the system. That 
behavior is in correspondence with the self-similarity of quasicrystalline structures, 
what allows scaling between the various hierarchic levels of the tilings. The resulting 
scaling exponents provide an insight into the properties of the macroscopic systems. 
Such a scaling approach was also used for the numerical determination of the electronic 
properties for the Penrose tiling, the Ammann-Kramer-Neri tiling, or octagonal tiling 
models 24,25!. However, all these calculations were performed for relatively small 
systems, making it difficult to determine the exact scaling behaviour. 

To obtain a better insight into the properties of quantum-mechanical wave 
functions and the quantum diffusion in quasicrystals, we study here models of d- 
dimensional quasicrystals with a separable Hamiltonian in a tight-binding approach 
[12] . This method is based on the Fibonacci sequences, which describes the weak and 
strong couplings of atoms in a quasiperiodic chain. The Fibonacci chain is the one- 
dimensional analogue of the Penrose tiling and the icosahedral quasicrystal 13 and, 
hence, it is related to quasicrystals with five-fold and ten-fold rotational symmetries. 
Our higher-dimensional tilings are then constructed as a direct product of these chains. 
Also the eigenstates of these tilings can be directly calculated by multiplying the wave 
functions and by adding/multiplying the energies of the quasiperiodic chains, which 
yields the hypercubic tiling or the labyrinth tiling, respectively. This approach allows 
us to study very large systems with up to 10 10 sites in three dimensions based on the 
solutions in one dimension. 

The quantum diffusion in the Fibonacci chain and the corresponding hypercubic 
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tiling has been studied quite intensively in the past [14H17] . Further, in our previous 
work we obtained more insight into the electronic eigenstates and the transport 
properties of the labyrinth tiling. For instance, we showed that the product structure 
of the labyrinth tiling leads to a connection of the multifractal properties of the wave 
functions in different dimensions [18,19]. In recent years also some new properties 
for one-dimensional quasiperiodic chains were found. These include the step-like 
behaviour of the wave-packet dynamics of metallic-mean chains reported by us |20j 
and the discovery of log-time-periodic oscillations of wave packets in the Fibonacci 
chain [21] . 

In this paper we report new results for the multifractal properties of wave 
functions in quasiperiodic systems and for the return probability of the wave packets. 
In particular, we study the participation numbers which provide an insight into the 
distribution of the probability density of electronic wave functions and allow us to 
distinguish extended and localized states. This can be related to the electronic 
transport properties, because e.g. Bloch functions indicate ballistic transport and 
localized wave functions are associated with insulators. To describe the electronic 
transport properties more quantitatively two quantities are usually studied: the mean 
square displacement and the temporal autocorrelation function. The scaling exponent 
f3 of the mean square displacement is connected to the electrical conductivity of a 
quasicrystal according to the generalized Drude equation [22 . In our recent work we 
showed that the scaling exponent (3 can be determined analytically for the labyrinth 
tiling in the regime of strong quasiperiodic modulation by a renormalization group 
(RG) approach [23]. In the present paper we focus on the temporal autocorrelation 
function, which is directly related to the properties of the (local) density of states. 
Further, it describes the decay of the wave packet which can be linked to the quantum 
diffusion of electrons in the system. We also show that for systems with absolute 
continuous energy spectra an improved lower bound for the exponent /3 can be 
obtained, which is not restricted to the systems studied here. 

The paper is structured in the following way. In Sec. [2] we briefly describe the 
construction rules for the d-dimensional separable tilings. We also introduce the 
participation numbers and the temporal autocorrelation function and give an overview 
of known results for these quantities in quasiperiodic systems. In Section [3] we briefly 
describe the RG approach, which is applied in the following two sections. This is 
followed by the presentation of new results for the properties of the wave functions 
in these systems in Sec. |4] We also apply the RG approach to show that the wave 
functions do not become localized even for very small coupling strengths. In Sec.[S]we 
study the wave-packet dynamics with a special focus on the return probability with 
numerical methods and with the RG approach. Further, the lower bounds for the 
scaling exponent j3 of the mean square displacement are discussed in Sec. [6] and the 
results are briefly summarized in Sec. [7] 

2. Construction, Electronic Structure and Wave-Packet Dynamics of 
Tiling Models 

2.1. Separable Tilings with Golden Mean 

The construction of the separable quasiperiodic tilings is based on the Fibonacci chain 
defined by the inflation rule V = {s — > w, w — > ws}. Starting with the symbol s we 
obtain after a iterations the ath order approximant C a of the Fibonacci chain. The 
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length f a of an approximant C a is given by the recursive rule f a = f a -i + f a -i with 
/o = /i = 1. Further, the ratio of the lengths of two successive approximants as well 
as the ratio of the numbers # of occurrence of the symbols w and s in an approximant 
approach the golden mean r for a — > oo. Thus, we obtain with the continued fraction 
representation r = [1] = [1,1,1,...] the relations 

lim = hm = t . (1) 

Our model describes an electron hopping from one vertex of a quasiperiodic 
chain to a neighboring one. The aperiodicity is given by the underlying quasiperiodic 
sequence of couplings, where the symbols w and s denote the weak and strong bonds 
of the chain. Solving the corresponding time-independent Schrddinger equation 

H = E l \¥) & ti-i^ti + ti,i+i*tf l = Ei *l (2) 
for the Fibonacci chain with zero on-site potentials, we obtain discrete energy values 
E % and wave functions \^ % ) = X^i" 1 10 represented in the orthogonal basis states 
\l) associated to a vertex I. If not stated differently we use free boundary conditions, 
which yields N a = f a + 1 eigenstates. The hopping strength t in the Schrodinger 
equation is given by the Fibonacci sequence C a with t s = s for a strong bond and 
t w = w for a weak bond (0 < w < s). For w = only isolated clusters of single 
vertices or two strongly coupled vertices occur, and for w = s the system is periodic. 

The (i-dimensional separable quasiperiodic tilings are then constructed from the 
product of d quasiperiodic chains which are perpendicular to each other. For this 
setup two special cases are known for which the systems become separable [12] : 

Hypercubic Tiling Hf d — This tiling corresponds to the usual Euclidean product 
of d linear quasiperiodic chains, i.e., only vertices connected by vertical and horizontal 
bonds interact as shown in Figure [TJ Such a tiling is referred to as a square tiling in 
two dimensions and in general as a hypercubic tiling [26,27]. In literature this model 
has been studied mainly for the Fibonacci sequence 15, 16,26-28-, and for reasons of 
comparison we will also present a few results of this model. 

Labyrinth Tiling — In the other case only coupling terms to neighbours along 
the diagonal bonds are considered (cf. Figure [lj, where the bond strengths of this tiling 
equal the product of the corresponding bond strengths of the one-dimensional chains. 
The resulting grid is denoted as labyrinth tiling. 

The Hamiltonians of these tilings are separable. Thus, the eigenstates of the 
tilings in d dimensions can be constructed from the eigenstates of d one-dimensional 



Wave Functions, Quantum Diffusion in Quasiperiodic Tilings 



5 



chains. In both cases the wave functions of the higher-dimensional tiling are 
constructed as the product of the one-dimensional wave functions, i.e., 

*;=*l;fci« (3) 

The superscripts s = (i,j,...,k) enumerate the eigenvalues E, and the subscript 
r = (l,m, . . . ,n) represents the coordinates of the vertices in the quasiperiodic tiling. 
However, the energies are treated in a different way for the two cases. In d dimensions 
the eigenenergies are given by 

E s = E l - J -- k = E u +E 2j + ... + E dk (4) 

for the hypercubic tiling and by 

E s = E iJ '-' k = E u E 2j ...E dk (5) 

for the labyrinth tiling. While for the hypercubic tiling all combinations of the indices 
i,j,...,k are allowed, for the labyrinth tiling only certain sets yield the correct solution 
due to symmetries of the eigenfunctions [18l[29] . 

2.2. Energy Spectra 

For the one-dimensional quasiperiodic systems the energy spectrum is known to be 
singular continuous [30H32] . The properties of the energy spectrum for the hypercubic 
tiling and the labyrinth tiling can be obtained from studying the scaling behaviour of 
the width of the energy bands [T21I27] . The numerical results for both tilings provided 
evidence for a transition from a singular continuous to an absolutely continuous 
energy spectrum upon increasing the coupling parameter w including an intermediate 
region where both parts coexist [12,27 ]. The transition parameters depend on the 
dimension and the underlying quasiperiodic structure. In particular, we like to point 
out the threshold values w 2 ^ ~ 0.6s in two dimensions and w 0.45s in three 
dimensions above which the spectrum of the hypercubic tiling and the labyrinth tiling 
becomes absolutely continuous [T21I2Z1I2II]- These transition parameters w t h seem 
to correspond to the threshold values wuos f° r which the density of states becomes 
gapless in two and three dimensions [33] . For w = the isolated clusters lead to a pure 
point spectrum and the periodic systems for w = s possess an absolutely continuous 
spectrum regardless of the dimension. 

Further, it was shown analytically that the energy spectrum of the square tiling 
is singular continuous in the regime of strong quasiperiodic modulation [31] . Another 
interesting property is that the energy spectrum of the Fibonacci chain splits into 
(sub)bands, which can be described by an RG approach proposed by Niu and Nori 
[34][35] . We use this property later to derive a relation between the structure and the 
properties of the electronic wave function. 

2.3. Participation Ratio and Multifractality of Wave Functions 

In Sec. 0] we study the multifractal properties of the wave functions in Fibonacci chains 
and the labyrinth tiling. This includes the computation of the participation ratio p, 
which reflects the fraction of the total number of sites for which the probability measure 
of the wave function | <I>^ | 2 is significantly different from zero. Hence, it provides 
information about the degree of localization of the wave function $ s [24]l36]. For 
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the discrete positions r the participation ratio is defined via the inverse participation 
number P($ s ) of a wave function <I> S and the number of sites V a in d dimensions as 

/ r ~\ — X\ 

P (Va) 



cx Vr 1 ■ (6) 



The expression (. . .) denotes the arithmetic average over all eigenstates in the spectrum 
or over a certain energy interval. In quasicrystalline systems the participation ratio 
p(&, V a ) of a state <£> was found to scale with a power law with respect to the number 
of sites V a , where the exponent 7( ( &) is connected to the nature of this eigenstate [24] . 
The wave function of a localized state is characterized by a scaling exponent 7 = 1, 
and 7 = corresponds to an extended state. Intermediate values of 7 (0 < 7 < 1) 
indicate fractal eigenstates, which are neither extended over the whole system nor 
completely localized at a certain position and show self-similar patterns [TU1EZ] - The 
participation ratios within a small energy interval usually fluctuate over a certain 
range, but the participation ratios and the inverse participation numbers as well as 
the corresponding scaling exponents do not show any significant trend over the whole 
energy spectrum for the systems considered here [19] , 

While the eigenstates of purely periodic systems can be characterized by a single 
scaling exponent D = D q , the wave functions of chaotic systems and quasicrystals arc 
usually multifractals [38-41]. The generalized dimensions D q describe the nature of 
the multifractality of the eigenstates, where D q is a monotonically decreasing function 
of q. The dimension D a is equal to the capacity dimension, which can be understood 
as the box counting dimension here, Di is equal to the information dimension, and 
D2 to the correlation dimension. 

The generalized dimensions can be determined either via the generalized inverse 
participation numbers or via a multifractal analysis. For the first method one 
investigates the scaling behaviour of the generalized inverse participation numbers 
averaged over all wave functions of a system with the linear system size N a , i.e., 

However, in Sec. I4.1l we present numerical results obtained via a multifractal analysis, 
which shows less numerical problems for negative parameters q in comparison to the 
generalized inverse participation numbers. The multifractal analysis is based on a 
standard box-counting algorithm, i.e., in d dimensions the system of linear system 
size N a is divided into B = N% / L d boxes of linear size L [4"2"H4l)] . For each box b we 
measure the probability to find the electron in this particular box, i.e., we study the 
spatial distribution of the wave functions. Hence, for an eigenstate $^ we determine 
the probability 

Mb ($ S ,i)= J2 l^l 2 (8) 
r£ box b 

and then compute the mass exponent 

e^o lne e^o his 

by a least squares fit of the quantity (\nP(q, $ s ,£)) versus lne considering box sizes 
from L = 10,..., N a /2 [USSE?]. The quantity e = L/N a denotes the ratio of the 
box size L and the system size N a . We only consider boxes with integer fractions 
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of N a /L 48J. In principle it is possible to improve the quality of the results by 
adapting the multifractal analysis to work also for non-integer fractions N a / L and 
by averaging over different initial positions of the boxes [451146) . However, we do not 
apply this improved method because it requires significantly more computing time. 
We also found that the simple partitioning scheme yields sufficiently accurate results 
for the determination of the mass exponents. The generalized dimensions D q are easily 
obtained from the mass exponents according to the relation D q = T q /(q — 1). 

2.4- Return Probability 

In the second part of the paper we investigate the dynamics of the quasiperiodic 
systems by means of the time evolution of a wave packet |T(i"o,t)) = 
^~J r6 £ T r (ro, t) |r), which is constructed from the solutions = ^e^ lE * of the 

time-dependent Schrodinger equation obtained by the separation approach. The 
wave packet is then represented as a superposition of these orthonormal eigenstates 
$r(*) where we assume that it is initially localized at the position ro of the 

quasiperiodic tiling, i.e., T r (ro,£ = 0) = S rro . Hence, with the completeness relation 
and normalized basis states the wave packet is defined by 

T r (r ,t) = J2K ^>- iESt . (10) 

s 

For a wave packet T r (ro,t) initially localized at the position ro the temporal 
autocorrelation function C(ro,t) corresponds to the integrated probability to be 
at the position ro up to time t > 0, i.e., 

C(r ,t) = l f \r ro (v Q ,t>)\ 2 dt> . (11) 

Here, we denote the integrand as return probability P(ro,t), i.e., 

F(r ,t) = |T ro (r ,t)| 2 . (12) 

The autocorrelation function for an infinite system is expected to decay with a 
power law C(r , t) oc t^ s ^ for t — > oo [49], where the scaling exponent S(r ) depends 
on the initial position ro of the wave packet. Further, it is known that the scaling 
exponent 6(ro) is equivalent to the correlation dimension D 1 ^ of the local density of 
states g(ro,E) of a system [4"9Tl52| . Due to this equivalence we can distinguish two 
limit cases and their associated spectral properties. For 6 — > the autocorrelation 
function and the return probability remain constant. A pure point spectrum implies 
£ — !• 0, but the opposite statement may not be true [33j[53]. In one dimension 6—^-1 
corresponds to the ballistic motion of an electron, where the electron moves freely on 
the chain without scattering [15]. In two and three dimensions 6—^-1 usually does 
not indicate ballistic transport. Further, an absolutely continuous spectrum implies 
S — > 1 but again the reverse is not necessarily true [53]. For quasiperiodic structures 
one often obtains intermediate values of S, i.e., < S < 1, which is related to the 
singular continuous energy spectrum [9HIUH5]. 

To provide a better insight into the dynamical properties and to compare the 
properties of different systems, we average the results over different initial positions of 
the wave packet, i.e., we investigate the scaling behaviour of C(t) = (C(r ,i)) cx t~ s . 
Previous studies mostly use non-averaged data, which do not adequately display the 
macroscopic transport properties [15|l29]. Due to restrictions in computing time we 
cannot consider all possible initial positions Tq. In particular, for the Fibonacci chain 
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we have chosen the 500 positions nearest to the center of the chain and in two and 
three dimensions we have randomly selected 100 vertices in a square of linear size 25 
and a cube of linear size 15 around the center of the tilings. 

Zhong and Mosseri already studied the temporal autocorrelation function and the 
return probability for the square and cubic Fibonacci tiling [15] . As for this model the 
higher-dimensional energy values are the sum of the one-dimensional energy values 
(cf. Eq. the time-evolution operator becomes separable. Therefore, the dynamics 
of the wave packets in higher dimensions are much better understood than for the 
labyrinth tiling. They were able to show that 

*id = S' M /d (13) 
and that the scaling behaviour of the temporal autocorrelation function is given by |15j 

c{t) « { i d s ld > i (14) 

Further, they found that for the limit case w = s one always obtains <5' ld = 1. With 
the same arguments, also the threshold values for the transition to an absolutely 
continuous energy spectrum can be directly obtained from the one-dimensional scaling 
exponents Sid according to d6id(w^ ) = 1. 



3. RG Approach 

In the regime of strong quasiperiodic modulation it is possible to relate the energy 
spectra and the wave functions of different approximants of the Fibonacci chain by an 
RG approach proposed by Niu and Nori [3H|54] ■ We use this approach to show in the 
following two sections that the wave functions do not become localized even for w — s- 0. 
In Sec. [5] we also apply the RG approach to find a relation between the quasiperiodic 
structure of the Fibonacci chain and its return probability. The general idea is that a 
given system can be transformed (renormalized) by the application of an RG step in 
such a way that we obtain a new system with a reduced number of degrees of freedom 
(i.e. sites/bonds). In particular, for the Fibonacci chain for w = the isolated sites 
(atoms) and biatomic clusters (molecules) coupled by a strong bond are visualized in 
Figure [5] This yields three highly degenerate energy levels: E = for the atomic 
sites and E — ±s for the bonding and antibonding states of the molecules [Mj- For 
non-zero coupling parameters w one finds a coupling between these isolated atoms and 
molecules, where the dominant contribution occurs between sites of the same type. 
Hence, for the Fibonacci chain C one can distinguish two possible RGs 34 : 



00000000000000000000<M>0000000<KWOOOO<KKXKW<>000 

C^^3 C^^3 O 0>=0 Ot=aO O C^^> 

O O^^^^^^— 



OO OOO OO O-O-O OOO O-O OOO O-O OOO O-O-O OO O-O-O OCK3 OO O-O-O O-O OOO O-O 



Figure 2. Atomic RG (top) and molecular RG (bottom) for the Fibonacci chain 
C: in the atomic RG all atomic sites survive an RG step, and in the molecular 
RG all molecular sites are replaced by a new site. The thick (thin) lines denote 
the strong (weak) bonds. 
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Atomic RG — The renormalized grid after one RG step consists of the atomic 
sites of the original chain as shown in Figure [5] These new sites are then connected 
by new strong and weak bonds. In particular, atoms separated by one molecule in the 
original chain become connected by a new strong bond and atoms separated by two 
molecules in the original chain get connected by a new weak bond. 

Molecular RG — In one RG step all molecular sites of the original chain are 
replaced by new atomic sites (cf. Figure [2]). The atoms in the renormalized grid are 
connected by new strong bonds between the sites of originally neighboring molecules 
and by new weak bonds between the sites for molecules that were separated by an 
atomic site. 

For both RG approaches we obtain a new Fibonacci chain, which is scaled in 
length and energy. The derivation of the length scalings c and the scaling factors z 
for the new bond strengths of an approximant under the application of one RG step 
is based on the Brillouin-Wigner perturbation theory. For the Fibonacci chain C one 
obtains for the [T31I54] 

• atomic RG a length scaling of c = r -3 and an energy scaling of z = w 2 / s 2 and 



• molecular RG a length scaling of c = t 2 and an energy scaling of z — w/2s. 

In other words, the energy spectrum £ a (s, w) of the Hamiltonian H (s, w) of the 
Fibonacci chain C a can be related to the spectrum of three Hamiltonians of smaller 
approximants [54] . With the scaling relations given above, it can be shown that the 
energy spectrum splits according to [M] 



This means that the central main band of the energy spectrum corresponds to the 
spectrum £ a -3(s,w) of the Hamiltonian H a _3(zs, zw) comprising the f a -3 atomic 
sites for w = scaled by a factor z. The left and right main bands at energies 
±s + z8 a ~2{s,w) correspond to the bonding and antibonding states of the spectra of 
the Hamiltonian H a -2(zs, zw) for the f a -2 molecular sites [55] . 

The recursive structure of the chains also has an influence on the properties of 
the individual wave functions. In particular, for the Fibonacci chain C a with periodic 
boundary conditions Piechon provided a quantitative result for the amplitudes of the 
wave functions for different approximants a [55] . Note that in this case the number of 
sites/eigenstates is N a = f a . For the energy levels E l (/ a _a < i + f a -2 < /a-i) in the 
central main band of the energy spectrum associated to the atomic sites, the relation 
is given by 



Likewise, one obtains the recursive relations for wave functions associated to the energy 
levels E % belonging to the left main band (1 < i < fa-2) and to the right main band 
(fa-i < i + fa-i < fa) of the energy spectrum: 



for the 




(15) 




(16) 




(17) 
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Figure 3. Scaling of the amplitudes of the wave functions \P| for the Fibonacci 
chain C a for the molecular RG (top) and the atomic RG (bottom) for w = 0.1 
and s = 1. 



The sites l + and l~ are those linked by a bond of strength s, which is replaced by 
the new site I' during one step of the molecular RG. These are approximate results 
because the atomic (molecular) RG approach ignores the probability of being on a 
molecular (atomic) site. Comparing this with the numerical results in Figure |31 we 
see that this nevertheless is a good approximation since the probability density on the 
other type of cluster is much lower in the regime of strong quasiperiodic modulation 
(w < s). 

Note that this RG approach is different from the RG approach used by Wang 
and Sanchez for the analysis of the electronic transport in Fibonacci chains and 
the corresponding hypercubic tilings pX 56 . With their renormalization rules they 
were able to derive an iterative procedure for the evaluation of Kubo-Greenwood 
formula at zero temperature. While this method allows one to efficiently compute the 
conductivity for large quasiperiodic systems, the RG approach applied in this paper 
allows us to directly connect the structure of the systems to certain properties of the 
wave functions and the transport properties. 

4. Multifractality of Wave Functions 

In this section we first present numerical results for the multifractal properties of 
the wave functions for the Fibonacci systems. We then show that the multifractality 
is directly related to the structure of the quasiperiodic systems for small coupling 
strengths w and that even for w — > the wave functions remain multifractals. 

4-1- Numerical Results for Multifractality of Wave Functions 

In Figure Q] we show numerical results for the scaling exponent 7 of the average 
participation ratio according to Eq. ([SJ for the Fibonacci chain and the labyrinth 
tilings. The data indicate that the scaling exponents for the labyrinth tiling in two and 
three dimensions approach the one-dimensional results for large system sizes [I8ll33j . 
We recently showed analytically that this property is fulfilled in the limit of infinite 
systems [19] . Additionally, for all systems the scaling exponents 7 approach a constant 
value for w — > 0. This limit behaviour can be related to the hierarchical structure of 
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Figure 4. Correlation dimension D2 of the inverse participation number Z2 and 
the scaling exponent 7 of the participation ratio p averaged over all wave functions 
of the systems for s = 1 . The largest considered systems are the Fibonacci chain 
C20 with 10 947 sites and the labyrinth tilings with 3 335 945 sites and £f % 
with 3 203 226 sites. The dotted lines in the insets correspond to the analytically 
predicted values D2 = 0.443 and 7 = 0.416. 



the chains and is further discussed in Sec. 14.21 

The multifractal properties of the wave functions are obtained by a multifractal 
analysis described in Sec. 12.31 The scaling behaviour of the quantity (In P(q, <£, L)) 
with the box size L is visualized for the Fibonacci chain in Figure [5] For most of 
the systems the results can be well described by a power law. Of course, for large 
negative values of q the accuracy deteriorates. A similar scaling behaviour can be 
observed for other values of w and for the higher-dimensional systems although we 
find some deviations from the power law for higher-dimensional labyrinth tilings and 
small coupling parameters w. The corresponding results for the generalized dimensions 
D q are shown in Figure |6] The results yield a box counting dimension of Dq = d, 
which corresponds to the spatial dimension d as expected. With increasing coupling 
parameter w the range of the scaling exponents D q narrows, i.e., the strength of the 
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Figure 5. Scaling behaviour of the function (In P(q, L)) used in the 
multifractal analysis with the ratio e for the Fibonacci chain C21 for ro = 0.1 
and s = 1. 
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Figure 6. Scaling exponent D q /d determined via a multifractal analysis for the 
Fibonacci chain C21 (solid lines) and the two-dimensional labyrinth tiling £jg 
(dashed lines) for s = 1. 



mult ifr act ality is stronger for smaller coupling parameters w. While in the limit w = s 
there is a single scaling exponent D q — d, we find that for w — > the distribution 
approaches a limit function, which can be related to the quasiperiodic structure. This 
is considered in more detail in the following section. 

Further, the results in Figure [5] indicate that the scaling exponents of different 
dimensions are related by D^ d — dD^ d . This is not surprising as we have shown 
analytically that for the golden-mean labyrinth this relation is approached with 
increasing system size |19) . 

4-2. Limit Behaviour of the Scaling Exponent 7 and the Generalized Dimensions D q 
for w — >• 

If the coupling strength w of the weak bond approaches zero, one can observe that the 
scaling exponent 7 of the participation ratio approaches a constant till it drops to zero 
for w = 0. This behaviour was already reported before |33j but no explanation for 
this phenomena was given. A similar behaviour can also be found for the generalized 
dimensions D q . This indicates that the wave functions do not become localized even for 
very small couplings w and that they remain multifractals for w — > regardless of the 
dimension. Thereby, this limit behaviour can be related to the hierarchical structure 
of the eigenstates of the quasiperiodic chains. In particular, for the Fibonacci chain C a 
we can use the recursive structure of the energy spectrum according to Eq. (|15j) and 
of the wave functions according to Eqs. (fT6]) and (fTT]) to derive the limit analytically. 

We start by deriving the participation ratios for an atomic state (f a -2 < i+fa-2 < 
fa-i) of the Fibonacci chain 
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Analogously, we obtain for the molecular states of the left main band of the energy 
spectrum (1 < i < fa-i) the expression 

-l 
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As these transformations do not directly depend on the energy we obtain the 
same expression also for the right main band of the energy spectrum. The average 
participation ratio is then directly calculated from these results according to 

Pa(fa)= ^EV^/aHr £ Pt™^ , fa) 

Ja <=1 hl i = U-2 + l 

1 fa 

— E Pa°\*\fa) (20a) 
i=/o-i+l 
/a-3 

J]pa-8(* i |/.-8) 



/a 

RG /a-3 1 



fa fa—c 



fa-2 1 



/a-2 



£pa-2(*\/a-2) 



(206) 

(20c) 



fa fa-2 ( 
= T _3 p a -3(/ a - 3 ) + T~ 2 p a _ 2 (/a-2) 

With the scaling behaviour of the average participation ratio p(f a ) /~ 7 according 
to Eq. (jSJ), we obtain the relation 

l =r -3(l+7) +T -2(l+7) . (21) 

This expression is independent of the coupling parameter w and yields a scaling 
exponent 7 ~ 0.416, which corresponds very well to the numerical results for w — > 
in Figure |4j 

With the same approach we can also derive an analytical expression for the 
generalized dimensions D q . Here, we start with the generalized inverse participation 
numbers for an atomic state according to Eq. ([7]), i.e., 

fa 
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Likewise, we obtain for the molecular states of the left main band of the energy 
spectrum 

z^m = Ei^i 29 - T ~ q E i*m 29 ( 23 «) 

1=1 Vemol 

fa-2 

= 2t-« £ 1*1' \ 2q = 2r-%, a - 2 (* 1 ) . (236) 
V=l 

The same expression is again valid for the right main band of the energy spectrum. 
This yields the recursive equation for the average generalized inverse participation 
numbers 

fa-2 fa-! 

z g , a = T J2 z T,a(^) + T E Z tr^ 1 ) 

<=1 »=/«-a+l 



1 



/a 



y- E ( 24 «) 

/a i=/a-l+l 
/a-3 



HG /a-3 1 _ r _ 3o/2 y 



/"a fa / _ | 



/a-2 



2 ^T i 7^ 2r ~ 9 E ^--2(^) (246) 

Ja Ja— 2 



i=l 
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With the scaling behaviour of the average generalized inverse participation numbers 
Zq(fa) <x fa T * according to Eq. (JT)), we obtain the relation 

I = r 3(r,-g/2-l) + 4r 2(T -g/2-l) _ (25) 

This expression is again independent of the coupling parameter w and yields e.g. 
for q — 2 the scaling exponent D 2 — — 0.443, which corresponds well to the 
numerical results for w — > in Figure [U While this analytical result shows an even 
better agreement with the numerical results for larger values of q, there are significant 
differences for smaller values of q as visualized in Figure [51 This is caused by the 
fact that the measures for negative parameters q are dominated by the boxes with 
the smallest probabilities. However, these probabilities are not considered in the RG 
approach because we neglect the probability on the other type of cluster. Nevertheless, 
the RG approach allows us to relate the quasiperiodic structure to the multifractal 
properties of the system for q > 2. 



5. Return Probability 

In order to obtain a deeper understanding of the connections between the transport 
properties and the quasiperiodic structure of a system, we study in this section the 
wave-packet dynamics for the different tilings. In the limit of weak coupling (w <C s) it 
is again possible to connect the structure of the systems and the wave-packet dynamics 
by an RG approach. 




Figure 7. Average temporal autocorrelation function C(t) for the Fibonacci 
chain C20 (left) and the associated labyrinth tilings C\§ (center) and (right) 
in two and three dimensions for s = 1. 



5.1. Scaling Behaviour of the Temporal Autocorrelation Function 

Some typical results for the temporal autocorrelation function C(t) introduced in Sec. 
12.41 are shown in Figure [7J for the Fibonacci chain and the associated labyrinth tilings 
in two and three dimensions. The corresponding scaling exponents S obtained by least 
squares fits are compiled in Figure [8j The power-law behaviour can be observed over 
several orders of magnitude in time, before C(t) approaches a constant due to finite 
size effects. While we obtain C(t) = 1 for a system with localized wave functions, 
we expect to get C(t — > oo) « 1/N a for a periodic chain of finite size N a due to 
the completely extended wave functions. A comparison with our numerical data for 
w = s = 1 in Figure [JJ shows a good agreement with the latter result. 

Further, the quantity C(t) decreases more strongly for higher values of the 
coupling parameter w in agreement with the more extended wave functions in these 
systems. However, for a periodic chain in one dimension one does not obtain the 
expected result S = 1 by a least squares fit due to sublogarithmic contributions [15 . 
Additionally, for small values of w we observe some oscillations. This corresponds to 
a step-like behaviour of C(r ,i), which can be related to the hierarchical structure of 
the chains 20J . 

Comparing the scaling behaviour of the average temporal autocorrelation function 
C(t) in FigureJJJwith respect to the dimension, we observe that for the one-dimensional 
chain the absolute value of the slope steadily increases with increasing coupling 
parameter w, while in two and three dimensions C(t) is nearly identical for w 2d > 0.6s 
and w 3d > 0.5s. These differences are also clearly visible for the corresponding scaling 
exponents 5 shown in Figured) which approach one when the coupling parameter w 
becomes larger than the threshold parameters lo?^ ~ 0.6s and « 0.45s for which 
the spectrum becomes absolutely continuous. It is expected that S — 1 for infinite 
systems if w > w?£. This limit behaviour of the scaling exponent 6 is indicated by 
the dashed lines in Figure [8j Similar results can be found for other metallic-mean 
systems [33] . Although 6—^1 indicates ballistic transport for the one-dimensional 
periodic chains [15] , the behaviour in two and three dimensions does not correspond to 
a transition towards ballistic spreading. The latter result is not always correctly used in 
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Figure 8. Scaling exponents <5 of the temporal autocorrelation function and 
the expected limit behaviour (dashed lines) for the infinite systems (top) as well 
as the scaling exponent 5' /d of the return probability P(t) (bottom) averaged 
over different initial positions of the wave packet for the Fibonacci chain and the 
labyrinth tilings in two and three dimensions. 

literature because e.g. Zhong and Mosseri wrongly state that there is a crossover from 
ballistic-like behaviour to a non-ballistic-like behaviour in the hypercubic tiling [15 . 
The results for the mean square displacement of the wave packet for the hypercubic 
or the labyrinth tiling clearly reveal anomalous transport for w < s. Only for w = s 
ballistic transport is observed (201125] . 

5.2. Scaling Behaviour of the Return Probability 

In order to obtain information about transport properties instead of the properties of 
the energy spectrum we can make use of the fact that already the integrand P(t) in 
Eq. (fTTj) shows a power-law behaviour according to P(t) cx t~ s . The integration itself 
is only used for smoothing the results because P(t) shows considerable oscillations (cf. 
e.g. Figs. [9l and [T0|. Naturally, one would expect that the scaling exponents of C(t) and 
P(t) are equal, i.e., 6 = 5'. However, the integration leads to two basic disadvantages. 
For 6' = 1 we obtain an additional logarithmic contribution C(t) cx \n(t)/t from the 




time t 

Figure 9. Moving average of the return probability P(t) over 10 data points 
(light colors) and the corresponding temporal autocorrelation function C(t) (dark 
colors) for the Fibonacci chain C20- Both quantities are averaged over 500 initial 
positions of the wave packet. 
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integration [15] . Hence, for the periodic one-dimensional systems and the commonly 
used system sizes one obtains exponents 6 significantly lower than 1 in agreement 
with our result 5 0.87. Further, in higher dimensions the scaling exponent 5' of 
P(t) becomes greater than 1 for large coupling parameters w (cf. Figure [5]). Zhong 
and Mosseri pointed out that in this case the integral in Eq. (ITT]) converges and yields 
a constant factor, which results in C(t) oc 1/t. Hence, the scaling exponent 5 = 1 
is obtained in the absolute continuous regime [15] . From Figure [8] we find 5' > 1 
for w > 0.55s in two dimensions and w > 0.43s in three dimensions in reasonable 
agreement with w?£ and w^. Since the integral converges rather slowly for coupling 
parameters w only slightly larger than Wth, we do not observe 6 = 1 for the finite 
systems considered here. However, for an infinite system we expect to obtain a clear 
transition towards 5 = 1 as indicated by the dashed lines in Figure [5] 

Therefore, we lose information about the transport properties and 6 = 1 is 
misleading because it often does not occur where it is expected (e.g. for periodic one- 
dimensional systems) but wrongly suggests a ballistic behaviour for w > u>th in two and 
three dimensions. Hence, for studying the transport properties it is better to extract 
the scaling exponent 6' via the calculation of P(t), or the temporal autocorrelation 
has to be redefined so that the integrand does not converge anymore. We choose the 
first option and use a moving average instead of the integration to smooth the results 
as shown in Figure [5] 

In Figure [8] the scaling exponent 6' of the return probability P(t) divided by 
the spatial dimension d is shown for the Fibonacci chain and the labyrinth tilings. 
The comparison of the scaling exponents for the temporal autocorrelation C{t) and 
its integrand P(t) shows that both results are quite close for the one-dimensional 
systems. The main difference occurs for w > 0.8s, where we obtain the expected value 
8' m 1 for w = s while the evaluation of C(t) results in 5 sa 0.87. Further, there 
are some differences for the Fibonacci chain, especially, for small values of w. This 
probably originates from the fact that we have not considered large enough times in 
the calculations of C(t). Actually, it is much more time-consuming to compute C(t) 
because the integration requires the consideration of much smaller time intervals. This 
is another advantage of calculating only the quantity P(t). 

Further, in two dimensions we observe that the scaling exponents S' 2d are about 
twice as large as the one-dimensional scaling exponents 5' ld . This is in good agreement 
with the results for the hypercubic tiling in Eq. (I13p . However, the results for the three- 
dimensional labyrinth tiling do not entirely fit into this scheme. While for w < 0.6s 
the expression 6' 3d /d is only slightly smaller than the exponent 6' ld , the behaviour 
changes completely for large values of w. In the latter case the scaling exponent 6' 3d 
becomes almost constant and approaches S' 3d (w = s) sa 1.75. The exact reason for 
this behaviour is unclear yet. Having a look at P(r , t) for different initial positions of 
the wave packet for w = s, one can observe that after an unusually smooth decrease 
of the function P(tq, t) according to a power law, there is an abrupt dip in the curves 
when finite size effects first occur (cf. Figure ITOl) . 

Since the square Fibonacci tiling and the two-dimensional labyrinth tiling are 
also identical for w = s, the scaling exponents 5 and 6' are identical. However, in 
three dimensions for w —> s the labyrinth tiling becomes a body centered cubic (bec) 
lattice while the cubic tiling approaches a simple cubic (sc) lattice. Therefore, these 
two grids have a different number of nearest neighbours independent of the coupling 
parameter w, i.e., eight in the labyrinth tiling and six in the cubic tiling. While the 
numerical results for the two-dimensional labyrinth tiling satisfy the relations in Eqs. 
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Figure 10. Scaling behaviour of the return probability P(ro, t) for the bee lattice 
with 240 3 /4 sites and the sc lattice with 240 3 sites as well as the analytical results 
of the return probability P(ro,t) for the bec lattice according to Eq. i'27ai . 



(|13p and (fT4|) , there are significant differences in the wave- packet dynamics in three 
dimensions, which are probably related to the different grid structure. Therefore, in 
Figure[TT]we plot for both systems the probability density of a three-dimensional wave 
packet along the y-z plane through the center of the x axis for different times t and the 
corresponding return probability P(ro,t) in Figure [TU] While the probability density 
of the wave packet for the sc lattice is distributed in a cube around the origin, for 
the bec lattice it spreads in the shape of an octahedron from the center with some 
prominent resonances along the x, y, and z directions. For the labyrinth tiling the wave 
packet already reached the boundaries for lni = 3 and was reflected due to the free 
boundary conditions. This leads to constructive interference at two boundaries and to 
destructive interference at the opposite boundaries because the initial position of the 
wave packet is chosen from one of the four possible center positions of the labyrinth for 
even chain lengths, which results in a slight asymmetry of the wave-packet dynamics. 

According to Figure [10] the return probability for the sc lattice strongly oscillates 
and decreases with a scaling exponent of 6' — 3, while for the bec lattice it decreases 
rather smoothly with a scaling exponent 6' « 1.75. The dip for \nt > 4 for the bec 
structure is clearly due to finite size effects and thus has not to be considered for the 
fitting of the scaling exponent. Although the center of the wave packet decreases more 
slowly for the bec structure than for the sc lattice, the overall wave packet spreads 
significantly faster for the bec lattice. This might be related to the higher number of 
nearest neighbours in the bec lattice. 

The data for the sc lattice are in good agreement with the results of Zhong and 
Mosseri (cf. Sec. I2.4|) . Their approach can also be used for the determination of the 
return probability P(t) — |T ro (i)| for the bec lattice. The approach makes use of the 
energy spectrum of the bec lattice in k space [57] 

£( k ) = 8coE > COf A cc . s ^ . (26) 

The probability amplitude at the initial position ro of the wave packet can be 
determined by the Fourier transform of this energy spectrum. Hence, we obtain 

T ro (r ,i) = fe~ iEt d^(E) = f e'^^dk (27a) 
J Jbz 
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Figure 11. Probability density along a plane through the center of the bee lattice 
with 240 3 /4 sites (left column) and the sc lattice with 240 3 sites (right column) 
at the times hit = 2 (top), licit = 2.5 (center), and lni = 3 (bottom). The 
corresponding values of P(ro,t) are shown in Figure flOl 
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We use the fact that the integration over the Brillouin zone (BZ) can be reduced 
to an integral over a cube from to 2ir due to the symmetries of the components 
of the k vector [57]. The function Jo denotes the Bessel function of the first kind. 
The integral over the Bessel function is evaluated for various time steps by numerical 
integration. The results agree exactly with the numerical data as shown in Figure 
1101 Unfortunately, it is not possible to evaluate this integral for arbitrary large times 
because the Bessel function strongly oscillates, which leads to numerical problems. 
While it is possible to derive with this approach an analytical expression for the long- 
time behaviour of the return probability for the sc lattice, we were not able to obtain 
a corresponding result for the bcc lattice. 



5.3. RG Approach for the Fibonacci Chain 

While we considered the limit case of the periodic systems in the previous section, 
we can also find a connection between the structure of the Fibonacci chain C and the 
transport properties in the regime of strong quasiperiodic modulation. In particular, 
we are able to derive an analytical expression for the scaling exponents 6' of the 
return probability P(t). This is based on the results of the RG approach, which yields 
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a scaling behaviour of the energy spectrum according to Eq. (jT5j) and of the wave 
functions according to Eqs. and (|T7|) . 

First, we determine the return probability P a (l ,t) for an approximant a for the 
two different kinds of initial positions ^o- For Iq being an atomic site we obtain 

natom/i +\ J -^ t l^atom/; jA|2 
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In the last step we make use of the fact that for large times t only the average behaviour 
of cos 2 (si) is relevant, which results in a factor of 1/2. The average return probability 
for an approximant a is then given by 
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With the scaling behaviour of the return probability P(t) oc t~~ 5 , we obtain the 
equation 

1 = t- 6 z- s ' +4t- 4 z~ s ' , (31) 

which approaches the numerical results for decreasing values of w as shown in Figure 
[T2l Note that Piechon [55] studied a similar quantity, which he also denoted as return 
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Figure 12. Comparison of the analytical and the numerical results of the scaling 
exponent 8' for the Fibonacci chain C20- 



probability. However, that quantity is related to the spectral measure of the density 
of states. In contrast, the return probability P(t) is related to the spectral measure 
of the local density of states for the initial position l of the wave packet. Only the 
latter quantity is related to the wave-packet dynamics. 

6. Lower Bound for Scaling Exponent of the Mean Square Displacement 

Another quantity often considered for the description of the electronic transport 
properties is the mean square displacement (also called width) 

i 

1 2 

d(r ,t) = 5>-r | 2 |T r (r ,t)| 2 (32) 
.rec 

of the wave packet. The spreading of the width d(ro,t) of the wave packet for an 
infinite system shows anomalous diffusion for t — >• oc according to d(rn,i) oc t^^ ra \ 
where the scaling exponent /3(rn) depends on the initial position ro of the wave 
packet [8,29,58-6(1. The electronic transport properties are governed by the wave- 
packet dynamics averaged over different initial positions, i.e., d(t) — (d(r ,t)) oc t@. 
Thereby, the scaling exponent ft is related to the conductivity a via the generalized 
Drude formula for zero-frequency conductivity a ~ e 2 g(Ep)ct^~ 1 [8122) . The quantity 
e denotes the elementary charge of an electron, g(Ey) the density of states at the 
Fermi level E-p, c a constant, and i sca a characteristic time beyond which propagation 
becomes diffusive due to scattering [SSI [HI]. Hence, ft — corresponds to the 
absence of diffusion, ft = 1/2 to classical diffusion, and ft = 1 to ballistic spreading. 
For quasiperiodic structures one often observes anomalous diffusion characterized by 
< ft < 1 [11,25,29 . For w <C s the width of the wave packet can be related by an 
RG approach to the structure of the Fibonacci chain [55] and the labyrinth tiling [23] . 

The computation of the scaling exponent ft is computationally very expensive 
especially in higher dimensions because it includes an average over all sites of the 
system. However, for the hypercubic tiling the scaling exponents are related according 
to ftdd — ftid due to the separability of the time-evolution operator. Numerical 
and analytical results suggest that this relation could be also valid for the labyrinth 
tilings [23,29,33 . However, in general one can make use of several lower bounds for 




Figure 13. Comparison of the scaling exponent ft with the lower bounds D^/d, 
D2/D2, and 5' / D2 for the Fibonacci chain (left) and the corresponding labyrinth 
tilings in two dimensions (center) and three dimensions (right). 



the scaling exponents /3, which have been proposed in literature. We shortly introduce 
two of them here and compare them to our numerical data [29,62 64 . Further, we 
present an improved lower bound for systems with absolute continuous energy spectra. 

Lower bound by Guarneri: As a rule of thumb, the wave-packet propagation is 
faster for smoother spectral measures. In particular, Guarneri proved for the scaling 
exponent j3 that in one dimension the lower bound j3 > D± holds 63; . This can be 
generalized for arbitrary dimensions to 

P > *f ■ (33) 

This expression quantitatively relates the diffusive properties to the information 
dimension of the spectral measure of the local density of states. A comparison with 
our numerical results in Figure Q2] shows that this inequality is clearly satisfied in one 
dimension. The spectral dimensions are computed by a multifractal analysis with 
respect to the energy. Therefore, we divide the energy range into B boxes of size AE 
and compute for a certain vertex ro of the tiling the function 

r(A£) = /W £ \<t>*A \k(AE) d 2 , (34) 

\b=l \£S£box& / / 

which shows a power-law behaviour with the generalized spectral dimension D% 
[381165]. 

With the inequality in Eq. (|33l) it is possible to show that an absolutely continuous 
energy spectrum in one dimension implies ballistic transport |66j and that in two 
dimensions this inequality requires a superdiffusive evolution of the wave packet 
(/3 > 5) in the presence of an absolutely continuous spectrum. However, in three 
dimensions it allows the occurrence of subdiffusive wave-packet spreading (/3 < |) 
although the energy spectrum is purely absolutely continuous, i.e., = 1 [55] ■ This is 
a very interesting aspect because investigations on real quasicrystals revealed that they 
are usually accompanied by absolutely continuous spectra [66] but the conductivity 
is well described by the generalized Drude formula for the case (3 < \ [81 1221 [5^1 15T] . 
Since D^ dd eventually becomes 1 for an absolutely continuous spectrum, this bound 
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strongly underestimates the actual scaling exponent /3 in two and three dimensions 
for large values of w as shown in Figure Q2J 

Lower bound by Ketzmerick et al.: Using the fact that the center of the wave 
packet is known [55] to decay according to t~ D * , Ketzmerick et al. showed that due 
to the normalization condition the spreading of a wave packet is described by an 
exponent j3 = Dt^/d in d dimensions as long as D 2 l < 1, i.e., for singular continuous 
spectra [51]. For wave packets which spread in a space of reduced dimension D 2 
instead of d dimensions they showed that the relation 
rC 

P = =f (35) 

holds for all systems which are characterized by a single scaling exponent 
/3 = f3q [Bl]. For the more general case of multiscaling in time (d q (t) = 
Sr^r l r ~ r °l 9 l^r(ro,i)| 2 oc t q/3 " ) they found that this yields a lower bound for 
positive moments, i.e., 



D 



^>-ff- ( 36 ) 

For negative moments one obtains a respective upper bound. 

Comparing this with our numerical results for q = 2 in Figure [T3] for the one- 
dimensional case, we find that the inequality in Eq. (|36[) is fulfilled and provides a good 
lower bound, which is significantly better than the result by Guarneri in Eq. (|33p (631 . 
The largest differences occur for w — > s, which are partly caused by numerical issues 
because the relation = 1 is satisfied for a periodic chain. The numerical results 



only fulfill the lower bound according to Eq. ([36]) because the wave-packet dynamics 
for the Fibonacci chain show multiscaling in time. This inequality also holds in two 
and three dimensions. However, in the regime of an absolutely continuous spectrum, 
the spectral dimensions fulfill — 1, and the inequality in Eq. (I36[) is no longer a 
good lower bound as shown in Figure H3l 

Addressing the latter case we would like to point out that the decay of the center 
of the wave packet is only described by the correlation dimensions of the local density 
of states for w < Wth because for w > Wth we obtain D% = 1 and the integral of the 
autocorrelation function converges [15) . However, in this regime the scaling exponent 
S' , which describes the decay of the center of the wave packet, can become larger than 
1 (cf. Figured]). Hence, for d-dimensional systems the bound in Eq. (|36|) should be 
replaced by 

5' 

P > 77- (37) 

in the absolutely continuous regime. According to Eq. ([13"]) for the hypercubic tiling 
this bound is as good as for the one-dimensional chain. This bound should hold also in 
general because is solely makes use of the normalization condition of the wave packet. 
We also checked in Figure [13] whether this relation is satisfied for the labyrinth tiling 
and found that this is a significantly better lower bound than that given in Eq. (|36[) 
for w > Wth- However, in three dimensions there are still large deviations compared 
to j3 for coupling parameters w — > s. This suggests that further studies are necessary 
in order to relate the width of the wave packet with the decay of its center. 

The benefit of good lower bounds for j3 is that the computation of the scaling 
exponents 5', D%, and D 2 is much less time consuming because for 5' it is not necessary 
to average over all grid positions r and the spatial and spectral dimensions, D 2 and 
Z?2 , are even time-independent. 
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7. Conclusion 

In this paper we presented the scaling exponents of the participation ratio and 
the inverse participation numbers for large approximants of the separable Fibonacci 
tilings. The results showed that the wave functions of the systems are multifractals and 
that the strength of the mult ifr act ality increases with decreasing coupling parameter 
w in agreement with the literature. By applying an RG approach we were able to show 
that the limit behaviour of the wave functions for w — > is related to the structure of 
the Fibonacci chain. Up to now this observation was only based on numerical results. 

Further, we presented numerical data for the quantum diffusion of electrons in 
these systems. Compared to previous results for the labyrinth tiling, we considered 
larger approximants and included an average over different initial positions of the wave 
packets, which allowed us to determine the scaling exponents with a higher precision. 
Our results also show that for the study of the quantum diffusion it is advisable 
to compute the return probability P(t) instead of the temporal autocorrelation 
function because for an absolutely continuous energy spectrum the integral over the 
return probability converges and no further information can be obtained. In higher- 
dimensional systems 5 — > 1 usually corresponds to the existence of an absolutely 
continuous part in the energy spectrum rather than to the occurrence of ballistic 
transport. 

In order to understand the diffusive properties of a system in detail it is also 
necessary to compute not only the width d(f) of the wave packet but also the return 
probability P(t). This can be easily understood from our results of the two tilings 
in three dimensions. While for the labyrinth tiling and the cubic Fibonacci tiling 
the scaling exponents /3 hardly differ for the two models [23], the scaling exponents 
5' are rather different for large coupling parameters w (cf. Figure I10|) . In Figure [TT] 
we have also seen that the wave packets spread in a quite different way for w — s. 
However, these differences are averaged out during the computations of the mean 
square displacement. 

For the Fibonacci sequence we were also able to derive an analytical expression 
of the scaling exponent 5' of the return probability P{t) in the regime of strong 
quasiperiodic modulation, i.e., we found a relation between the quasiperiodic structure 
of the chain and the emerging electronic transport properties. Further, we show that 
the exponent 5' can be used to define a better lower bound for the scaling exponent (3 
of the width of the wave packet for systems with absolutely continuous energy spectra. 
This result is not restricted to the quasiperiodic systems considered here. 

References 

[1] Penrose R 1974 Bull. Inst. Math. Appl. 10 266-271 

[2] Shechtman D, Blech I, Gratias D and Cahn J W 1984 Phys. Rev. Lett. 53 1951-1953 

[3] Fujiwara T 1999 Physical Properties of Quasicrystals ed Stadnik Z M (Springer, Berlin 

Heidelberg New York) pp 169-207 
[4] Senechal M 2006 Notices of the American Mathematical Society 53 886—887 
[5] Levine D and Steinhardt P J 1984 Phys. Rev. Lett. 53 2477-2480 

[6] Dubois J M and Weinland P 1993 Coating materials for metal alloys and metals and method 

European Patent EP 0356287 Al and US Patent 5204191 
[7] Stadnik Z M 1999 Physical Properties of Quasicrystals (Berlin Heidelberg New York: Springer) 
[8] Roche S, Trambly de Laissardiere G and Mayou D 1997 J. Math. Phys. 38 1794-1822 
[9] Damanik D 2000 Directions in Mathematical Quasicrystals vol 13 ed Baake M and Moody R V 
(American Mathematical Society, Providence) pp 277-304 
[10] Kohmoto M, Sutherland B and Tang C 1987 Phys. Rev. B 35 1020-1033 



Wave Functions, Quantum Diffusion in Quasiperiodic Tilings 



25 



[11] Sire C and Mosseri R 1989 J. Phys. France 50 3447-3461 

[12] Sire C 1989 Europhys. Lett. 10 483-488 

[13] Socolar J E S and Steinhardt P J 1986 Phys. Rev. B 34 617-647 

[14] Abe S and Hiramoto H 1987 Phys. Rev. A 36 5349-5352 

[15] Zhong J and Mosseri R 1995 J. Phys. : Cond. Matt. 7 8383-8404 

[16] Lifshitz R 2002 J. Alloys Compd. 342 186-190 

[17] Sanchez V and Wang C 2004 Phys. Rev. B 70(14) 144207 

[18] Thiem S and Schreiber M 2010 J. Phys. : Conf. Ser. 226 012029 

[19] Thiem S and Schreiber M 2011 Eur. Phys. J. B 84 415-421 

[20] Thiem S, Schreiber M and Grimm U 2009 Phys. Rev. B 80 214203 

[21] Lifshitz R and Even-Dar Mandel S 2011 Phil. Mag. 91 2792-2800 

[22] Schulz-Baldes H 1997 Phys. Rev. Lett. 78 2176-2179 

[23] Thiem S and Schreiber M 2012 Phys. Rev. B 85(22) 224205 

[24] Grimm U and Schreiber M 2003 Quasicrystals: Structure and Physical Properties ed Trebin 

H R (Wiley- VCH, Berlin) pp 210-235 

[25] Passaro B, Sire C and Bcnza V G 1992 Phys. Rev. B 46 13751-13755 

[26] Ilan R, Liberty E, Even-Dar Mandel S and Lifshitz R 2004 Ferroelectrics 305 15-19 

[27] Even-Dar Mandel S and Lifshitz R 2008 Philos. Mag. 88 2261-2273 

[28] Ashraff J A, Luck J M and Stinchcombe R B 1990 Phys. Rev. B 41 4314-4329 

[29] Yuan H Q, Grimm U, Repetowicz P and Schreiber M 2000 Phys. Rev. B 62 15569-15578 

[30] Damanik D and Gorodctski A 2011 Comm. Math. Phys. 205 221-277 

[31] Damanik D and Gorodctski A 2009 Electron. Res. Announc. Math. Sci. 16 23-29 

[32] Siito A 1989 J. Stat. Phys. 56 525-531 

[33] Ccrovski V Z, Schreiber M and Grimm U 2005 Phys. Rev. B 72 054203 

[34] Niu Q and Nori F 1986 Phys. Rev. Lett. 57 2057-2060 

[35] Niu Q and Nori F 1990 Phys. Rev. B 42 10329-10341 

[36] Wegner F 1980 Z. Phys. B: Condens. Matter 36 209-214 

[37] Repetowicz P, Grimm U and Schreiber M 1998 Phys. Rev. B 58 13482-13490 

[38] Brandes T, Huckestein B and Schweitzer L 1996 Ann. Phys. 508 633-651 

[39] Guarneri I and Mantica G 1994 Phys. Rev. Lett. 73 3379-3382 

[40] Huckestein B and Klesse R 1997 Phys. Rev. B 55 R7303-R7306 

[41] Salejda W 1995 Int. J. Mod. Phys. B 9 1453-1474 

[42] Bottgcr H and Kasner G 1990 physica status solidi (b ) 162 489-495 

[43] Castcllani C and Pcliti L 1986 J. Phys. A 19 L429 

[44] Chhabra A and Jensen R V 1989 Phys. Rev. Lett. 62 1327-1330 

[45] Grussbach H and Schreiber M 1993 Chemical Physics 177 733-741 

[46] Schreiber M and Grussbach H 1991 Phys. Rev. Lett. 67 607-610 

[47] Vasquez L J, Rodriguez A and Romer R A 2008 Phys. Rev. B 78 195106 

[48] Rodriguez A, Vasquez L J and Romer R A 2008 Phys. Rev. B 78 195107 

[49] Ketzmerick R, Petschel G and Geisel T 1992 Phys. Rev. Lett. 69 695-698 

[50] Barbaroux J M, Germinct F and Tcheremchantsev S 2001 Duke Math. J. 110 161-193 

[51] Gucrin C A and Holschneider M 1997 J. Stat. Phys. 86 707-720 

[52] Holschneider M 1994 Comm. Math. Phys. 160 457-473 

[53] del Rio R, Jitomirskaya S, Last Y and Simon B 1995 Phys. Rev. Lett. 75 117-119 

[54] Piechon F, Benakli M and Jagannathan A 1995 Phys. Rev. Lett. 74 5248-5251 

[55] Piechon F 1996 Phys. Rev. Lett. 76 4372-4375 

[56] Sanchez V, Perez L A, Ovicdo Roa R and Wang C 2001 Phys. Rev. B 64(17) 174205 

[57] Jclitto R J 1969 J. Phys. Chem. Solids 30 609-626 

[58] Schulz-Baldes H and Bellissard J 1998 Rev. Math. Phys. 10 1-46 

[59] Poon S J 1992 Adv. Phys. 41 303-363 

[60] Huckestein B and Schweitzer L 1994 Phys. Rev. Lett. 72 713-716 

[61] Mayou D 2000 Phys. Rev. Lett. 85 1290-1293 

[62] Barbaroux J M, Combes J M and Montcho R 1997 J. Math. Anal. Appl. 213 698-722 

[63] Guarneri I 1993 Europhys. Lett. 21 729-733 

[64] Ketzmerick R, Krusc K, Kraut S and Geisel T 1997 Phys. Rev. Lett. 79 1959-1963 

[65] Halsey T C, Jensen M H, Kadanoff L P, Procaccia I and Shraiman B I 1986 Phys. Rev. A 33 
1141-1151 

[66] Bellissard J and Schulz-Baldes H 2000 J. Stat. Phys. 99(1) 587-594 

[67] Mayou D, Berger C, Cyrot-Lackmann F, Klein T and Lanco P 1993 Phys. Rev. Lett. 70 3915- 
3918 

[68] Guarneri I 1996 J. Math. Phys. 37 5195-5206 



